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Abstract 

A high-degree spherical harmonic series is used to compute the radial 
deformation of the Earth by oceanic tidal loading. By exploiting “fast” 
numerical transforms, this approach is found to be much more efficient — 
but no less accurate — than the traditional Green’s function approach. The 
method is used to derive an atlas of load tide maps for 10 constituents of the 
NSWC ocean tide model. 
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Introduction 


Tidal deformations of the Earth affect a wide range of geophysical and geode- 
tic measurements {Baker, 1984; Lambeck, 1988). This is particularly true for 
the modern, high-precision techniques of space geodesy: satellite and lunar 
laser ranging, very long baseline interferometry, and satellite altimetry are 
all seriously affected by the deformation in their range measurements and/or 
their inferred baselines to other stations. Our purpose here is to discuss one 
particular type of deformation — radial deformation due to ocean loading — 
and to present an atlas of maps displaying this deformation for some of the 
major tidal constituents. 

Our primary motivation in producing such an atlas is to provide an ad- 
ditional correction to satellite altimeter measurements. The range of the 
loading tide can easily reach 10 cm, and this is far too large an error source 
to ignore — compare, for example, the stringent requirements of the TOPEX 
error budget as listed by Bom et al. (1984). 

In satellite altimetry, only the geocentric tide ( g can be directly observed. 
This tide consists of a body tide &, a land-relative ocean tide ( „ , and a load 
tide Ci'- 

Cg = Cb + Co + Cl- 

The body tide is relatively well understood and requires for its accurate eval- 
uation only the astronomical, tidal-generating potential (given by Cartwright 
& Edden (1973)) and the Love number h 2 . (h 3 is also required for some mi- 
nor lunar components of the generating potential.) The ocean tide £„ is more 
problematic. There are a number of numerical models of Co available; none 
is completely satisfactory, and, in fact, the improvement of our knowledge 
of Co is a major goal of altimetry. Probably the best current model of 
for the deep-ocean regions is that computed at the Naval Surface Weapons 
Center (NSWC) by Schwiderski (1980, 1983); this model tide is provided on 
the GEOSAT geophysical data records as a standard tide correction for the 
altimetry. The load tide Ci is a measure of the deformation of the solid Earth 
under the weight of the ocean tide. Ci as presented below corresponds to (is, 
in fact, derived from) the NSWC ocean tide. 

Direct tidal analysis of satellite altimetry is already promising to improve 
our knowledge of the ocean tide ( Cartwright & Ray , 1989). In tidal analysis 
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of altimetry, it is usually more direct to derive an “altimetric tide” £ a , which 
is neither a geocentric tide nor a land-relative ocean tide but is given by 

Ga = Gg ~ Gb = Go Cl- 

There are several points to make here. If the goal is simply to provide 
accurate tidal corrections to the altimetry so that other, more subtle oceano- 
graphic signals may be studied, £a is actually what is required. If only £ a is 
available, as is the case for current GEOSAT data processing, then Gi must 
be computed to form £ a . But in addition to this use, knowledge of £t in its 
own right has several applications. Among these are: (a) to compare the 
satellite tide £ a with tide gauges which are direct measurements of (b) to 
calibrate an altimeter by using local observations of sea level; (c) to compute 
£ 0 from £ a in order to provide corrections for Earth-tide measurements of 
gravity, tilt, and strain ( Harrison , 1985); and (d) as a specific application, 
to estimate the loading tide contribution to the C„(GEOSAT)— £ 0 (NSWC) 
difference maps of Cartwright & Ray (1989). 

The usual procedure for computing tidal loading effects for Earth-tide 
observations is to integrate the ocean tide over the surface of the globe using 
a numerical Green’s function, as in the classic studies by Longman (1963) 
and Farrell (1972). This integration is moderately time consuming, but since 
there are normally only a few Earth-tide stations, the time is of little concern. 
For satellite altimetry, however, this is no longer the case, since the load must 
be evaluated at very nearly the same resolution as the ocean tide, which is 
1° X 1° for the current NSWC model. Efficient methods are therefore required. 
Thus, our secondary motivation in this work has been to experiment with the 
applicability of replacing the Green’s function methodology with one based 
on a high-degree-and-order spherical harmonic expansion. This has allowed 
us to exploit recent advances in algorithm design for “fast” transforms. 

The atlas of computed load-tide charts is presented in Appendix C. A 
previous suite of global loading maps was published by Parke & Hendershott 
(1980) and Parke (1982). Their maps are derived from their own ocean-tide 
model, which has somewhat fewer constituents and lower spatial resolution 
than the NSWC model. Quite recently, Francis (1989) has presented results 
similar to our own by a different method. 
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Formalism 


In this section, let the tide £ be a complex number whose real and imaginary 
parts are the in-phase and quadrature parts of a particular tidal constituent. 
Similarly, let U be a complex potential. 

We will now briefly recapitulate the formalism for loading tide calcula- 
tions. The theory is standard (e.g., Lambeck , 1988), but it is useful to have 
it summarized. The starting point, following Munk & MacDonald (1960), is 
to write the radial deformation C i as a combination of loading Love numbers 
that are defined by 

G = £*n WJg, (1) 

n 

where U' n is the n-th degree spherical harmonic component of the potential 
U' caused by the mass of the ocean tide For a spherical Earth of radius 
® e > Co, U' is given by 

U\0 , <f>) = Gp w a\ J f dtt 

with p w the density of seawater, r is the distance between (0 , <f>) and ( 9 ', <f>') 
and, as usual, may be expressed as 


1 _ 1 
r 2a e siny>/2 


— fl p n(cosiJ>) 


c n=0 


and, by the addition theorem, as 


( 2 ) 


; = f £ E 

r a e n =Om=-n + 1 

where the Y nm are normalized, complex spherical harmonics. The tide 
may be expressed in spherical harmonics: 


«M)=£ £ 

n=0 m=— n 


where the coefficients are 


( 3 ) 

( 4 ) 
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which allows the potential to be expressed as 

W(t,+) = i*Gp.a. £ £ (5) 

n=0m=-n Zn ' L 

The loading tide may now be expressed in two alternative ways. Combining 
(1) and (5) gives 

oo n out 

9 r 1 a nm^nm(^) <f>) (6) 

Pe n=0 ro=— n 2n + 1 

with p e the mean density of the Earth. Alternatively, combining (1) and (2) 


gives 

= p^\ J 1 U0'A')GU’)dO! 

(7) 

where 

oo 



W) = TT 2 kn Pn( COS V>) 

m « n - 0 

(8) 


and M e is the mass of the Earth. 

The two approaches (6) and (7) may be loosely described as the “wavenum- 
ber domain” and “space domain” approaches, respectively. In the past, 
both approaches have been used for computing loading corrections for var- 
ious Earth-tide measurements, where, depending on the measurement type, 
the “operator” 3h|,/2n-f 1 or the Green’s function G(i>) is modified as re- 
quired. Workers who have used the wavenumber domain for loading calcula- 
tions include Pertsev (1966), Groten & Brennecke (1973), and Goad (1980). 
Nonetheless, more success has been achieved using (7) rather than (6), the 
primary reasons being efficiency and complications related to slow conver- 
gence of the series. The series (3) converges relatively slowly, which requires 
the maximum degree, say AT, of the series to be large and thus the number 
of coefficients o nm , evaluated by the surface integrals (4), to be quite high. 
On the other hand, after the function G(if . >) has been tabulated, (7) requires 
the evaluation of only one surface integral at each Earth-tide station. In our 
present application, however, we require (i to be evaluated over the entire 
globe (or at least, over all the oceans). Below, we describe some algorithms 
that allow the summations in (3) and the discrete version of (4) to be carried 
out efficiently for large N (for 1° resolution, N = 180 is appropriate). The 
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result is that approach (6) will be found to be much more economical than 
(7) for our present application. 

Furthermore, computing the loading effects of radial deformation is much 
more amenable to the use of spherical harmonics than computing any of 
the other quantities normally required in Earth-tide studies. The operator 
Zh , n /2n -\- 1 in (6) ensures that the series for Ci converges faster than the series 
for ( 0 . In contrast, the operator to compute loading corrections for gravity 

n + — (n + jX . . 

2n + l ^ ' 

which causes the series for gravity to converge no faster than the series for 
Co, even allowing for the fact that the Love number k' n decays as ~ 1/n (see, 
however, the modification of (9) due to Merriam (1980)). The operators 
for tidal tilt and strain are worse yet, causing convergence more slowly than 
that for Co • We may therefore expect that the spherical harmonic approach 
to quantities other than radial deformation may be less successful, although 
we have done no further investigation of the matter. 

2.1 Spherical harmonics and fast transforms 

When progressing from two-dimensional Fourier series on a rectangle to 
spherical harmonics on a sphere, it must seem a natural extension to attempt 
to incorporate the fast Fourier transform (FFT) algorithm. The matter has 
been studied by a number of authors since Ricardi & Burrows (1972), but 
the extension is only partially successful. Useful summaries are by Swarz- 
trauber (1979) and Colombo (1981). As they note, it is clear that, during 
evaluation of the coefficients a nm of (4), the integration over <f> can be han- 
dled by an FFT. There exists, however, no fast Legendre transform for the 
integration over 0, so at this point a number of methods have been proposed. 
Some authors (e.g., Brown , 1985) have simply performed standard numerical 
quadratures for the Legendre transforms; others (e.g., Dilts , 1985) have ex- 
pressed the Legendre functions as Fourier series so that a second FFT could 
be invoked; other approaches have been used as well. As Colombo (1981) has 
noted, the most efficient methods require 0( N 3 log IV) arithmetic operations 
compared to 0(JV 4 ) operations for straightforward quadrature of (4). (For 
comparison, Fourier analysis on an N x N square requires but 0(N 2 log N ) 


5 



operations.) 

The procedures we have used in our loading calculations are based on 
the work of Swarztrauber (1979) (see also Browning et al ., 1989). He uses an 
auxiliary coefficient function Z™(6) which allows the harmonic coefficients to 
be computed by 

2 N N 

«nm = E 

j=l i= 1 

The functions Z™ for spherical harmonic analysis thus correspond to the 
Legendre functions P™ for spherical harmonic synthesis. The reader may 
consult the references for further details. 

3 Accuracy Assessments 

We divide our discussion on the accuracy of the loading tide calculations 
into two sections. The first concerns the computational method: spherical 
harmonics vs. Green’s functions. The second section concerns a more general 
assessment independent of the computational method. 

3.1 Comparisons over the Indian Ocean 

This section reports on some comparisons of the spherical harmonic and 
Green’s function approaches to ocean loading calculations of radial deforma- 
tion. This testing is important for several reasons. Firstly, even though we 
are using a high-degree spherical harmonic series (N = 180), it is important 
to establish that this is indeed high enough. After all, Farrell (1972) used 
N = 10000 for evaluating his Green’s function using the series (8). Secondly, 
it is natural to expect certain manifestations of the Gibbs’ phenomenon — at 
least in the synthesis of by (3) — particularly near coastal areas where the 

ocean tide jumps discontinuously from a possibly large value to zero over 
land. The effect of this on the radial deformation must be monitored. 

In Figures 1, 2, and 3 are shown the results of three loading tide com- 
putations for the Ma constituent in the Indian Ocean area, derived from the 
(global) Ma model of Schwiderski (1983). Figure 1 was computed using the 
Green’s function tabulated by Farrell (1972) for the Gutenberg-Bullen Earth 
model. (Some further details of this are given in Appendix A.) Figures 2 and 
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tide, using 
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tide, using 



3 show the results of the spherical harmonic series approach, with N = 36 
and N = 180, respectively, using the Gutenberg-Bullen Love numbers h' n 
also tabulated by Farrell (1972). 

From Figure 2 it is apparent that the resolution provided at N = 36 is 
too coarse, even though the large-scale structure is probably adequate. The 
agreement between Figures 1 and 3 is quite remarkable, and establishes the 
validity of the spherical harmonic approach. Although there is some “wig- 
gliness” in the contours over Africa in Figure 3, they are of no consequence; 
most importantly, there is no serious indication of any Gibbs-type ringing, 
even near Madagascar where the NSWC tide is over lm in amplitude. 

Figure 3 (as well as Figure 2) was actually part of a global calculation — it 
in fact is a portion of the global Figures C5-6 of Appendix C. To produce 
this global figure required almost exactly 10 times less computer-time than 
did Figure 1, which was computed only for the area shown in the figure. 

3.2 General accuracy assessments 

The previous section established the accuracy of the spherical harmonic cal- 
culation of radial deformation by comparing it to the standard Green’s func- 
tion calculation. There are several other issues that affect accuracy that are 
independent of the spherical harmonic/Green’s function choice. 

Ocean tide model. The most obvious issue concerns the accuracy of the 
ocean tide ( 0 — in our case, the NSWC numerical model. The NSWC M a 
model is advertised as accurate to better than 5 cm anywhere in the open 
ocean ( Schwideraki , 1983). A recent comparison of the Mj model with 66 
open ocean (island or bottom pressure) measurements shows a standard de- 
viation of 4.2 cm ( Cartwright & Ray, 1989), which would indicate the 5-cm 
limit is often exceeded. In certain locations the model is known to be inac- 
curate due to lack of real data. Nevertheless, as noted below, a comparison 
between our load tide maps (Appendix C) and the Parke-Hendershott (1980; 
Fig. 13) load tide maps shows good agreement, almost always within 1 cm — 
and a large fraction of that must be due to errors in the Parke-Hendershott 
ocean tide. This fact, and the 7% rule of thumb (i.e., Ci « — ' 0.07£ o , see 
below), would indicate that a reasonable estimate of the loading error due to 
errors in the ocean tide is 0.5 cm (in the open ocean). 
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In coastal and shelf areas, of course, the NSWC model is generally poor 
because the 1° resolution is too coarse to model such complex tides. There- 
fore, our corresponding load tide maps are likewise suspect in these areas. 
For accurate loading computations in such areas, the standard recourse (e.g., 
Baker , 1980) is to complement the open-ocean tide model with a local, high- 
resolution model. This proviso is particularly relevant to VLBI measurements 
from antennae located in coastal areas like, for example, the Gulf of Alaska. 

Mass conservation. In theory any model of the ocean tide should conserve 
mass. In the spherical harmonic decomposition (3) we should then have 
o 0 o = 0. For a variety of reasons, however, most numerical models of the tides 
violate this constraint, and some controversy has developed on how best to 
allow for this. The degree of non-conservation for the NSWC model may be 
seen in the list of spherical harmonic coefficients tabulated in Appendix B. 

Farrell (1972b) has stressed the need for forcing mass conservation in the 
ocean model to improve the residuals of gravity tides, although (in Farrell , 
1973) he reports that doing so by (a) removing a (complex) slab of water 
from the model or (b) forcing aoo = 0 gave similar results for the perturbed 
potential, an issue later clarified by Agnew (1983). We observed that sup- 
plementing Farrell’s set of Love numbers (which had no n = 0 term) with 
h' 0 = —0.134 from Longman (1963) results simply in a (bias) shift of (i by 
0.34 — 0.23i mm. This, as well as the 7% rule applied to the aoo terms of 
Appendix B, leads us to believe that mass non- conservation in the NSWC 
model is of little concern for our deformation calculations. 

Earth models. If we compare the load tide results for M a as given in the 
atlas with similar calculations that use the Love numbers of Pertsev & Ivanov 
(tabulated in Melchior, 1983) as well as of Zschau (1978), we find maximum 
vector differences of 1.39 and 1.42 mm, respectively, of which about 30% is 
due to the simple n = 0 bias term. We conclude that the choice of the gross, 
radially symmetric, frequency-independent Earth model is of little concern. 

Of more concern, however, may be lateral heterogeneities of the Earth 
that are not included in the Earth models. Effects of local geology or to- 
pography, for example, are known to play a large role — as much as 50% or 
more ( Farrell , 1979) — in Earth-tide observations of tilt and strain, less so 
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for gravity. The effect was first noticed when scientists compared nearby ob- 
servations that were unexpectedly discordant. Radial deformation is not (as 
of yet) directly measured at the required precision, so the effects of lateral 
heterogeneities on it are not completely evident. Finite element modeling, 
however, by Gong et al. (1975) and Beaumont (1978) indicate that effects 
larger than 10% are possible. In light of this, our maps should be inter- 
preted as the “homogeneous” load tide, following the terminology of Berger 
& Beaumont (1976). 

We should also point out that, even though it is easily incorporated into 
the spherical-harmonic methodology, we have not allowed for any frequency 
dependence in the Love numbers. In particular, the nearly- diurnal free wob- 
ble, due to the liquid outer core, affects the diurnal band near our Ki con- 
stituent by increasing h' n for the (n,m) = (2,1) harmonic only (Wahr & 
Sasao , 1981). (Resonance effects on the ocean model can be assumed to be 
incorporated automatically by the model’s observational constraints.) From 
Wahr & Sasao’s Table 5, the effect is seen to be less than 1% of the body 
tide for Ki, which, at this level, we may neglect. 

4 Discussion 

Appendix C presents an atlas of maps showing the amplitudes and Greenwich 
phases for the load tides corresponding to all but one constituent of the 
NSWC ocean model — namely the semidiurnals Kj, S 2 , M 2 , N 2 , the diurnals 
Ki, Pi, Oi, Qi, and the long-period tides Mf and Mm- All were computed 
using the Gutenberg-Bullen Love numbers (with no frequency dependence) as 
tabulated by Farrell (1972). (The NSWC model also includes the semiannual 
Ssa tide. We have not included load-tide maps for Ssa, since this component 
is influenced by large thermal and atmospheric effects that do not load the 
lithosphere.) 

The general features of the maps are in good agreement with the 6 charts 
previously published by Parke & Hendershott (1980) and Parke (1982), al- 
though, as one would expect, many of the details are different. In general, for 
the semidiurnal tides, the Parke- Hendershott amplitudes are slightly larger 
than our amplitudes, the largest discrepancy being nearly 10 mm in the M 2 
tide for the high-amplitude region in the Indian Ocean. This arises because 


12 



the Parke-Hendershott ocean tide there is nearly 10 cm larger than the NSWC 
tide. 

For the important M 3 constituent, we have converted the amplitude 
and phase plots to in-phase and quadrature components (i.e., HcosG and 
if sin G). This is for direct comparison to the Co(GEOSAT)— £„(NSWC) 
(color) difference maps of Cartwright & Ray (1989) and allows us to estimate 
the contribution of the load tide to this difference. These plots are shown 
in Figures 4 and 5. Comparing these to the Cartwright/Ray maps shows 
little similarity. Furthermore, the small amplitudes of the load tide — 4 cm 
or less — would not markedly change the (GEOSAT— NSWC) discrepancy ar- 
eas. One region of correlation is the large quadrature anomaly off the coast of 
northeastern Brazil; but the 3-cm load tide still leaves a (GEOSAT— NSWC) 
difference of nearly 15 cm. 

4.1 Coherence of ocean and load tides 

Various tidal authorities (e.g. Schwiderski, 1983) have used in the past the 
following rule-of-thumb for estimating the load tide Cl from the land-relative 
ocean tide Co- 

Cl = -eCo 

with e typically around 0.07 — thus, the “7% rule” referred to earlier in this 
work. 

Such a rule is obviously but a rough approximation. To see the extent 
of the approximation, we have computed the coherence between the NSWC 
ocean tides and the load tides of Appendix C for a number of constituents. 
The coherence is given by (with C complex) 

T , KC6)! 2 

7 (Co’Co) (CfCi) 

where ( ) denotes averaging over the global ocean. The results are summa- 
rized in Table 1, along with the least-squares estimates of e. e is shown as 
a real number; allowing it to be complex gives an insignificant imaginary 
component. All the e estimates for the short-period tides are less than 7%. 
Performing the calculations only over the deep (> 1km) ocean improves 7 s 
marginally and increases e at most ten percent. 
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Table 1: Coherence and scaling coefficients. 


Tide 

7 2 

6 

M a 

0.84 

0.047 

s a 

0.86 

0.046 

Ki 

0.83 

0.057 

Oi 

0.87 

0.062 

M f 

0.94 

0.078 

M m 

0.86 

0.066 


In Table 1, the e- values for the short-period tides agree almost exactly 
with similar values computed by Parke (1982) for the Par ke-Hender shot t 
ocean tide, except for Ki, which is here somewhat smaller. These values of e 
are not comparable with the larger values of Accad & Pekeris (1978), since 
they include (in their variable K) effects of ocean self-attraction, which do 
not concern us here. 

Finally, to show even more plainly the limitations of this “7% rule,” we 
have used our Appendix C maps to compute the vector error in the ap- 
proximation for e = 0.0667 (the value used by Schwiderski (1983)). This is 
shown as Figure 6. ( Francis (1989) has recently computed a similar map, 
with nearly identical results.) The errors are clearly largest near coastlines; 
in the Gulf of Alaska, for example, the error in near 100%. Figure 7 shows 
the errors when e = 0.047. The large errors near coastlines are somewhat 
reduced, at the expense of increased errors in the open ocean where they 
tend to mimic the amplitudes of Q shown in Figure C5, indicative of e being 
too small. The error in the mid-Pacific approaches 15 mm — nearly 50% of £/. 
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Appendix A Integrating Green’s Functions 

We discuss in this appendix a few details concerning the calculation of loading 
tides with a (numerical) Green’s function. 

The loading tide, given by Eq. (7), is discretized as follows: 

Cl(0> <t>) = Pv>a\ J Co(0', <f>') G(ip) dO! 

= (ai) 

since the ocean tide is a constant over each 1° x 1° block ft;. The Green’s 
function G(ip) for radial deformation has an (integrable) 1/ip singularity as 
ip — ► 0 (Farrell, 1972), which complicates evaluation of the integral in (Al). 
However, when ip is sufficiently large, the integral can be approximated using 
the centroid rule 

/ f a Q(1>) dO’ = G(ipi) A ft; = Q(ipi) sin 0[ A0[ A ft 

where ip, is the distance from ($, <p) to the center of 0, and Aft; is the solid 
angle subtended by ft;. Numerical tests indicate this approximation yields 
an error typically no more than 5% when ip > 1° and no more than 1% when 
iP > 5 °. 

For the blocks ft; in the immediate neighborhood of the computation 
point (6, (p), a more accurate integration rule must be used. We have used for 

this region an adaptive quadrature algorithm due to van Dooren & de Ridder 
(1976). This algorithm does not require the evaluation of the integrand on 
the boundary of the region, and it can therefore handle a singularity there. 
Thus for the block ft; which contains the point (0, <p) where ip = 0, we split 
the block into two regions with (0, <p) on the boundary and integrate the two 
regions separately. 

Calculations of the contribution of individual regions to the summation 
in (Al) as a function of ip shows a quite complex and unpredictable behavior, 
due to the complicated behavior of the tides. It is therefore essential to carry 
out the integration over the entire globe even though G(ip) decays rapidly 
away from the origin. 
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Appendix B Tide Coefficients 


As a by-product of our spherical harmonic analyses of the NSWC ocean tide 
models, we list here the first few harmonic coefficients for each constituent. 
They may be of some interest in their own right. To interpret them, how- 
ever, requires knowledge of the normalization, and it seems nearly everyone 
differs in this regard. As noted elsewhere, we have used a procedure due to 
Swarztrauber, who expresses a spherical harmonic series as (all quantities 
real) 

TV TV n 

\ ^2 a n0 P°(cos 6) + '£2 cos m<j> - b nm sin m<j>)P ™ (cos 6). (Bl) 

n=0 n=lm=l 


The normalization is set such that P ^(m) * s normalized when it (alone) is 
integrated from fi = — 1 to 1. Thus, 



dO = 7 r. 


A particular tidal constituent £ = H cos(cr£ — G ) with speed a and Green- 
wich phase G is decomposed into in-phase and quadrature components and 

Ca ; 

£i = H cos G (2 = H sin G. 


These components are then expanded in spherical harmonics: 


N 


Ci = H (°nm COS m<f> - b nm sin m<f>) P™( cos 0) 


n=0 m=0 


(2 = ^2^2 (^m cos m< t > - dnm sin TTl<f>) P™( COS 6) 

n=0 m=0 

where the prime on the second summation indicates that the first term is 
multiplied by 1/2, in agreement with (Bl). The coefficients ®nm i b nm , Cjim j 
and d nm are listed in the following tables up to degree and order 6; units are 
millimeters. 

To convert our tabulated coefficients to the more conventional notation 
used in satellite geodesy (e.g. Lambeck , 1988), use the following formulae: 

Dnm COS £ nm = f ^ d nm )N n 

Dnm S ^ n £ nm = | (Cnm ~F b nm )N n 


22 


where the normalization factor is 


a rm 2n + l (n - m)! l 1/2 
n [ 2 (n + m)! 
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Coefficients of the NSWC tide K2. Coefficients of the NSWC tide S2. Coefficients of the NSWC tide M2. 
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Coefficients of the NSWC tide Oi. Coefficients of the NSWC tide Qi. Coefficients of the NSWC tide Mf. 
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Appendix C Atlas of Load Tide Maps 

On the following pages is an atlas of co-amplitude and co-phase maps for 
10 of the load tides corresponding to particular constituents of the NSWC 
ocean tide model. These constituents are (in order of increasing period) 
K 2 , S 2 , M 2 , N 2 , Ki, Pi, 0 1? Qi, Mf, M m . As noted above, the S sa con- 
stituent is not included in this atlas, since a substantial fraction of this tide 
is non-gravitational and non-loading. The maps were produced using the 
Gutenberg-Bullen Love numbers as tabulated by Farrell (1972). All maps 
are on the Miller cylindrical projection. 
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